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Abstract  -  This  paper  investigates  the  optimization 
of  both  single  and  full  polarization  radar  transmission 
waveforms  to  maximize  target  identification  discrimina¬ 
tion.  This  theory  is  applied  to  the  discrimination  of  the 
T-72  and  Ml  battle  tanks  based  upon  simulated  target 
frequency  response  data.  Significant  performance  im¬ 
provement  in  identification  is  obtained  using  an  opti¬ 
mized  transmission  waveform  over  that  of  a  standard 
chirped  pulse. 

1  INTRODUCTION 

A  number  of  researchers  [1,  2,  3,  4,  5,  6,  7]  have  com 
sidered  the  use  of  sophisticated  pulse  shaping  techniques  in 
order  to  maximize  the  radar  energy  reflected  off  of  a  non¬ 
point  target.  In  particular.  Grieve,  Guerci,  Pillai,  Oh,  and 
Youla,  [1,  2,  3]  have  developed  a  general  theory  of  opti¬ 
mized  pulse  shaping  that  maximizes  the  target  SINR,  in¬ 
cluding  the  effects  of  both  generic  colored  noise  and  col¬ 
ored  signal-dependent  clutter.  In  addition,  Guerci  and  Pillai 
[8,  9]  developed  the  theory  of  optimized  pulse  shaping  for 
single-channel  target  identification  discrimination  via  the 
use  of  techniques  that  are  similar  to  that  used  for  detection. 
This  paper  extends  this  target  discrimination  analyses  in 
permitting  multiple-chaimels,  colored  noise,  and  non-zero 
colored  clutter. 

The  present  analysis  applies  the  theory  of  optimized 
pulse  shaping  for  target  identification  discrimination  us¬ 
ing  two  simulated  surface  targets:  the  T-72  and  Ml  main 
battle  tanks.  SAIC-Champaign  [10]  generated  the  full- 
polarization  VHF-band  radar  signatures  for  a  single  eleva¬ 
tion  angle  of  15°  and  the  full  spectrum  0°  -  360°  of  aspect 
angles  relative  to  the  sensor.  These  VHF-band  data  were 


generated  using  the  Fast  Illinois  Solver  Code  (FISC)  that 
applies  a  method-of-moments  technique  to  provide  high  fi¬ 
delity  at  relatively  low  radar  frequencies.  The  specific  VHF- 
band  data  generated  by  SAIC-Champaign  cover  frequencies 
between  225-375  MHz  at  an  aspect  interval  of  2°. 

2  OPTIMIZED  SINGLE-POLARIZATION  TAR¬ 
GET  IDENTIFICATION 

The  derivation  begins  with  the  result  that  the  maximiza¬ 
tion  of  the  probability  of  correct  classification  between  two 
target  classes  a  and  is  equivalent  [12,  13]  to  the  maxi¬ 
mization  of  the  square  of  the  Mahalanobis  distance 

=  (ya  -  y/3)^R"Hya  -  yp)  (i) 

between  the  two  target  echoes.  Here,  =  qaf  and 
=  q^f  are  real-valued  vectors  of  length  2N  —  1  giving 
the  temporal  samples  of  the  echoes  from  targets  a  and  /?,  re¬ 
spectively.  The  real- valued  vector  f  —  [/o,  /i, . . . , 
gives  the  temporal  samples  of  the  transmission  pulse.  The 
real-valued  matrices  and  q^  are  the  convolution  impulse 
responses  for  targets  a  and  ^5,  respectively,  having  the  form 
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is  the  temporal  autocorrelation  of  the  noise  plus  clutter,  with 
matrix  coefficients 

re  =  -^f  {Gn(a;)  +  dw.  (4) 

Thus,  Tf  can  be  expressed  in  the  form 

7?2  ^  (5) 


with  the  matrix  defined  by 


n  =  (q^  -  q/3)^R  ^{qa  -  q/?)-  (6) 


6)  Let  ffc+i  ^  Fk+iiu;)  and  go  back  to  Step  2  with  k 
replaced  by  A:  +  1,  and  repeat  until  Ck  is  sufficiently 
small.  Then  the  optimized  transmission  vector  is 

f  =  lim  ffc.  (9) 

k—*OQ 

Figure  1  gives  the  improvement  in  the  square  of  the  Ma- 
halanobis  distance  squared  between  the  T-72  and  the  Ml  at 
VHF-band  resulting  from  the  use  of  the  optimized  transmis¬ 
sion  pulse  over  that  of  a  standard  chirped  pulse.  This  figure 
shows  two  values  of  the  CNR:  0  and  10.  The  improvement 
in  the  square  of  the  Mahalanobis  distance  degrades  as  the 
CNR  level  is  increased,,  as  occurs  with  the  SINR  improve¬ 
ment  in  the  detection  problem  [2]. 

For  the  case  of  aspect  uncertainty,  it  is  necessaiy  to  com¬ 
pute  the  expected  value  of  the  square  of  the  Mahalanobis 
distance,  i.e., 

^  =  Jd9  =  jde  (lO) 


For  the  case  of  zero  clutter  Gc(ct>)  =  0,  the  minimax  theo¬ 
rem  implies  that  the  maximization  of  rp-  is  obtained  when 
the  transmission  pulse  vector  f  is  equal  to  the  eigenvector 
of  O  corresponding  to  the  largest  eigenvalue. 

For  the  case  of  non-zero  clutter  Gc{oj)  ^  0,  the  au¬ 
tocorrelation  matrix  R  depends  upon  the  power  spectrum 
|jP(a;)p  of  the  transmission  pulse  vector  f  via  Eqs.  (3)  and 
(4),  so  that  an  iterative  procedure  similar  to  that  used  for 
optimized  target  detection  [2]  must  be  applied,  as  described 
below: 

1)  For  the  initialization  A:  =  0,  begin  with  any  real  causal 
temporal  vector  fo  of  duration  fo  and  energy  Eq, 

2)  Let  ffc  Fk(u))  and  find  the  corresponding  temporal 
autocorrelation  matrix  Rfc  using  Eqs.  (3)  and  (4). 

3)  Compute  the  ftk  matrix  using  Eq.  (6)  in  terms  of  the 
autocorrelation  matrix  Rfc  and  the  target  impulse  re¬ 
sponse  matrix  q. 

4)  Find  the  largest  eigenvalue  and  corresponding 
normalized  eigenvector  of  the  f^fc  matrix, 

5)  Define  the  error  at  stage  k  by 


and  invoke  the  same  update  rule  that  is  applied  in  Pillai 
and  Guerci  [2] 


with  the  density  function  ^(0)  characterizing  the  a  prior 
likelihood  of  the  target  aspect  0.  The  matrix  O(^)  now  in¬ 
cludes  aspect  dependence,  i.e., 

nw  =  (w„(0)  -  W;5(e))"R-'(WaW  -  w^(0)).  (11) 
Inserting  Eq.  (1 1)  into  Eq,  (10)  implies  that 

^  (12) 

can  be  expressed  in  terms  of 

n  =  y  d9^(e)^i{e).  (13) 

Thus,  optimization  of  the  transmission  waveform  to  maxi¬ 
mize  identification  performance  under  conditions  of  aspect 
uncertainty  involves  the  computation  of  the  weighted  aver¬ 
age  of  ft  (6)  matrices  with  respect  to  aspect.  Furthermore, 
the  iterative  procedure  described  above  for  the  case  of  non¬ 
zero  clutter  is  modified  only  by  the  replacement  of  the  O 
matrix  by  its  weighted  average  ft. 

3  OPTIMIZED  FULL-POLARIZATION  TARGET 
IDENTIFICATION 

This  section  describes  the  theory  of  optimal  waveform 
transmission  and  reception  in  order  to  maximize  the  Ma¬ 
halanobis  distance  between  two  target  echoes  for  the  case 
of  a  single  full-polarization  waveform,  i.e.,  one  contain¬ 
ing  both  horizontal  and  vertical  components.  Consider  the 
2A’-length  real-valued  transmission  signal  vector  and  cor¬ 
responding  frequency  response  vector 
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with  f/i,  ft;,  Ffi  and  each  containing  N  temporal  sam¬ 
ples.  The  subscripts  h  and  v  denote  the  horizontal  and  ver¬ 
tical  channels,  respectively.  This  transmit  vector  is  further 
constrained  to  have  finite  energy  Fq-  This  energy  constraint 
corresponds  to  the  case  in  which  the  sum  of  the  transmis¬ 
sion  energies  in  both  the  horizontal  and  vertical  channels 
are  fixed,  so  that  a  single  power  supply  supports  both  trans¬ 
mission  chaimels.  The  2N  x  2N  target  impulse  response 
matrix  and  corresponding  frequency  response  matrix  have 


the  form 

<lvv  J 

V  Wrfc 

) .  (15) 

The  target  echo  vector  has  the  form 

S  = 

(16) 

The  full-polarization  matrix 
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is  the  temporal  autocorrelation  of  the  noise  plus  clutter,  with 
the  2  X  2  sub-matrix  coefficients 


ri  —  —J  {Gn(c4;)  Gp{co)}€^^  (Lj.  (18) 

The  matrices  Gn(a;)  and  Gir(a;)  are  the  fiill-polarization 
spectral  densities  corresponding  to  the  noise  and  the  clutter, 
respectively.  The  total  clutter  power  spectral  density  has  the 
form 

Gf{uj)  =  G;,;,(a;)|F;,(a;)|2  +  G;,,(u;)F;.(a;)F;(a;)  (19) 

+  Gv/i(u;)Ft;(t*;)F;J'(£*;)Gvv(ci;)|Ft;(a;)p  >  0  (20) 

The  optimization  of  the  transmission  vector  f  in  order  to 
maximize  the  square  of  the  fiill-polarization  Mahalanobis 
distance  gives 

rp  =  (21) 

with  the  matrix  ft  defined  by 

=  {qa  -  -  q^}-  (22) 

For  the  case  of  zero  clutter  Gc(c«^)  =  0,  the  minimax  the¬ 
orem  implies  that  the  maximization  of  is  obtained  when 
the  transmission  pulse  vector  f  is  equal  to  the  eigenveaor 
of  ft  corresponding  to  the  largest  eigenvalue.  For  the  case 
of  non-zero  clutter  Gc(a?)  /  0,  the  autocorrelation  matrix 
R  depends  upon  the  full-polarization  power  spectrum  of  the 
transmission  pulse  vector  f  via  Eq.  (20),  so  that  the  iterative 


procedure  described  for  single-polarization  above  must  be 
applied. 

Figure  2  gives  the  full-polarization  waveforms  optimized 
to  maximize  the  Mahalanobis  distance  between  the  T-72 
and  the  Ml,  as  a  function  of  the  relative  aspect  angle  for  the 
case  of  white  noise  and  zero  clutter.  The  optimized  wave¬ 
form  typically  focuses  the  majority  of  its  energy  into  a  nar¬ 
row  frequency  band  corresponding  to  the  maximum  target 
response  at  that  aspect  angle.  Rgure  3  demonstrates  that 
the  optimized  full-polarization  waveform  gives  an  improve¬ 
ment  of  1-5  dB  in  the  Mahalanobis  distance  over  that  ob¬ 
tained  from  the  transmission  of  a  full-polarization  chirped 
waveform. 

The  analysis  described  above  for  the  case  of  aspect  cer¬ 
tainty  can  be  extended  to  the  case  of  aspect  uncertainty  in  a 
manner  similar  to  that  performed  for  the  single-polarization 
case.  The  resulting  theory  requires  a  weighted-average  with 
respect  to  relative  aspect  angle  be  performed  on  the  auto¬ 
correlation  kernel  matrix  ft.  This  averaging  of  the  full- 
polarization  kernel  matrices  yields  a  smoothing  of  the  Ma¬ 
halanobis  curves,  as  was  obtained  in  the  single-polarization 
case. 

4  CONCLUSION 

This  study  investigates  the  optimization  of  a  single  trans¬ 
mission  pulse  shape  and  the  receiver  impulse  response  in 
order  to  maximize  the  probability  of  correct  identification 
between  two  target  classes.  The  optimization  of  the  trans¬ 
mission  pulse  shaping  in  order  to  maximize  target  identifi¬ 
cation  performance  that  was  developed  by  Guerci  and  Pillai 
[9]  is  extended  to  include  multiple  channels,  colored  noise, 
and  non-zero  colored  clutter.  These  extensions  [11]  for  the 
identification  problem  are  developed  via  a  maximization  of 
the  Mahalanobis  distance,  and  thus  the  probability  of  cor- 
rea  classification,  between  the  echoes  of  two  target  classes. 

This  study  applies  this  theory  [9]  and  extensions  of  op¬ 
timized  transmission  pulse  shaping  in  order  to  investigate 
the  maximization  of  the  probability  of  correct  identifica¬ 
tion.  Algorithmic  implementation  for  the  simulated  T-72 
and  Ml  frequency  response  data  at  both  single  and  multiple 
polarizations  of  the  VHF  frequency  band  reveals  significant 
improvements  in  the  Mahalanobis  distance  of  using  a  single 
optimized  waveform  over  that  of  a  standard  chirped  pulse. 
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mission  pulse  shape  optimized  for  T-72  versus  Ml 
identification  discrimination.  The  two  curves  corre¬ 
sponding  to  CNR=  1  and  CNR=  10  are  plotted  as  a 
function  of  aspect  for  0°  —  360®. 
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Opt  waveforms  In  FroQuoncy  (H):  T72  -  M1^^.  »  0,  WR  *  0 


u. 


375Lt-„f_: — . . : — : - ilJ 

0  $0  180  270  360 

Aspect  Anglo  (dog) 

Figure  2.  This  figure  gives  the  full-polarization  wave¬ 
forms  optimized  to  maximize  the  Mahalanobis  dis¬ 
tance  between  the  T-72  and  the  Ml,  as  a  function  of 
the  relative  aspect  angle  for  the  case  of  white  noise 
and  zero  clutter.  The  optimized  waveform  typically 
focuses  the  majority  of  its  energy  into  a  narrow  fre¬ 
quency  band  corresponding  to  the  maximum  target 
response  at  that  aspect  angle. 


Mahalanobis  Distance  Comparison;  T72  -  Ml.o^  =  0,  CNR  =  0 
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Figure  3.  This  figure  demonstrates  that  the  optimized 
full-polarization  waveform  gives  an  improvement  of  1- 
5  dB  in  the  Mahalanobis  distance  over  that  obtained 
from  the  transmission  of  a  full-polarization  chirped 
waveform. 
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